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Abstract 

Hamiltonian systems of hydrodynamic type occur in a wide range of applications including 
fluid dynamics, the Whitham averaging procedure and the theory of Frobenius manifolds. 
In 1 + 1 dimensions, the requirement of the integrability of such systems by the generalised 



hodograph transform implies that integrable Hamiltonians depend on a certain number of 



arbitrary functions of two variables. 

On the contrary, in 2 + I dimensions the requirement of the integrability by the method 
of hydrodynamic reductions, which is a natural analogue of the generalised hodograph trans- 
form in higher dimensions, leads to finite-dimensional moduli spaces of integrable Hamilto- 
nians. In this paper we classify integrable two-component Hamiltonian systems of hydrody- 
namic type for all existing classes of differential- geometric Poisson brackets in 2D, establish- 
ing a parametrisation of integrable Hamiltonians via elliptic/hypergeometric functions. Our 
approach is based on the Godunov-type representation of Hamiltonian systems, and utilises 
a novel construction of Godunov's systems in terms of generalised hypergeometric functions. 
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1 Introduction 

A 1 + 1 dimensional quasilinear system, 

u t = A(u)u x , 

is said to be Hamiltonian if it can be represented in the form Uj = Ph u where h(u) is the 
Hamiltonian density, h u is its gradient, u = (u , . . . , u n ), and P = P tJ is a Hamiltonian operator 
of differential-geometric type, 

It was demonstrated in |10| that the coefficients g %i define a flat contravariant metric, and b l - k can 

be expressed in terms of its Levi-Civita connection via 6^ = — g^T^. This implies the existence 
of flat coordinates where the operator P assumes a constant-coefficient form, P tJ = e l S i 3 4-. In 
the same coordinates, the corresponding matrix A is the Hessian matrix of h. It was conjectured 
by Novikov that a combination of the Hamiltonian property with the diagonalisability of the 
matrix A implies the integrability. This conjecture was proved by Tsarev in [37J who established 
the linearisability of such systems by the generalised hodograph transform. Hamiltonian systems 
arise in a wide range of applications in fluid dynamics and the Whitham averaging theory. The 
underlying differential-geometric structure make them the main building blocks of the theory 
of Frobenius manifolds |12| . Let us note that the requirement of the diagonalisability of the 
matrix A imposes a complicated system of third order PDEs for the Hamiltonian density h 
whose general solution is known to be parametrised by n(n — l)/2 arbitrary functions of two 
variables. Thus, the moduli space of integrable Hamiltonians in 1 + 1 dimensions is essentially 
infinite-dimensional. We refer to |10| [38j [11] for further discussion. 

In this paper we provide a complete classification of integrable two-component 2 + 1 dimen- 
sional Hamiltonian systems of hydrodynamic type. Our approach is based on the method of 
hydrodynamic reductions which, in some sense, can be viewed as a natural extension of the 
generalised hodograph transform to higher dimensions. Our results demonstrate that the moduli 
spaces of integrable Hamiltonians in 2 + 1 dimensions are finite-dimensional: the integrabil- 
ity conditions for the corresponding Hamiltonian densities constitute over-determined involutive 
system of finite type. Modulo natural equivalence groups, this leads to finite lists of integrable 
Hamiltonians parametrised by elliptic/hypergeometric functions. 

Thus, we consider 2 + 1 dimensional quasilinear systems 

u t = A(u)u x + B(u)u y (1) 

which can be represented in the Hamiltonian form, 

u t = Ph u . (2) 

Here h(u) is the Hamiltonian density, h u is its gradient, and P = P %1 is a Hamiltonian operator 
of differential-geometric type, 

P ij = ^'(u)^ + #(u)«* + g*(u)-jL + 6«(u)tj{. 

Operators of this form are generated by a pair of metrics g %3 , g lJ , see [101 l25l I26| for the general 
theory and classification results. The main difference from the one-dimensional situation is that, 



although both metrics g lJ and g 10 must necessarily be flat, they can no longer be reduced to 
constant coefficient form simultaneously: there exist obstruction tensors. The two-component 
situation is understood completely. One has only three types of Hamiltonian operators: the first 
two of them can be reduced to constant-coefficient forms, 



d/dx 
d/dy 



(3) 



P=[ ° d Jff ), (4) 



and 

P -f . , 

d/dx d/dy J ' 

while the third one is essentially non-constant, 

p= /2vw\d_ + /0v\d_ + /v x v y \ _ 

\ w J dx \ v 2w J dy \ w x w y J ' 

Here v,w are components of the vector u, and h u = (h v , h w Y. Hamiltonian systems generated 
by the operators ([H]) - (0) take the form 

v t = {hv)x, Wt = (h w ) y , (6) 

Vt = (h w ) x , w t = (h v ) x + (h w ) y , (7) 

and 

v t = (2vh v + wh w - h) x + (vh w ) y , w t = (wh v ) x + (2wh w + vh v - h) y , (8) 

respectively. Our main result is a description of Hamiltonian densities h(v, w) for which the 
corresponding systems © - (|8|) are integrable by the method of hydrodynamic reductions as 
proposed in |13| . Let us point out that the integrability conditions for the general class of two- 
component systems ([1]) were derived in |14| in the coordinates where the first matrix A is diagonal. 
It was demonstrated in |27l [29] that, again in special coordinates, the matrices A and B can be 
parametrised by hypergeometric functions. However, it is not clear how to isolate Hamiltonian 
cases within these general descriptions. In this paper we adopt a straightforward approach and 
derive the integrability conditions by directly applying the method of hydrodynamic reductions 
to the systems ([6]) - ([8]). For Hamiltonian systems of type ([6]) this was done in [13j. Thus, in the 
Appendix we illustrate the method of hydrodynamic reductions by calculating the integrability 
conditions for equations of type 0. Applied to equations © - ©, this method results in 
involutive systems of fourth order PDEs for the Hamiltonian densities h(v, w) which are quite 
complicated in general. It was observed that these systems simplify considerably under the 
Legendre transformation h(v, w) — > H(V, W) defined as 

V = h v , W = h w , H = vh v + wh w — h, Hy = v, Hw = w. 

In the new variables, equations (JEJ) - (EJ) take the so-called Godunov form, 

(H v ) t = V x , (H w ) t = W y , (9) 

(H v ) t = W x , {H w ) t = V X + W y , (10) 

and 

(H v )t = (VH v + H) x + (WH v ) y , (H w ) t = (VH w ) x + (WH w + H) y , (11) 



respectively. All our classification results will be formulated in terms of the Legendre-transformed 
Hamiltonian densities H{V, W). We recall that a system ([TJ is said to be in Godunov form |2T] 
if, for appropriate potentials F a (u), a = 0, 1,2, it possesses a conservative representation 

{F ,i) t + (Fx,i) x + {F 2 ,i) y = 0, 

F a ,i = dF a /du % . In particular, in the case (fTTj) one has Fq = —H,F\ = VH,F2 = WH. This 
representation will be utilised in Sect. 5 to provide a parametrisation of the generic integrable 
potential of type (jlip by generalised hypergeometric functions. 

Let us first review the known results for systems of type ([9]) . Here the integrability conditions 
take the form 

Hy W H V vvv = ZFtvvvHvvWi 

HvwHyvvw = ZHwvHvww, 

HywHvvww = HyywHvww + HyyvHwww, 

HywHvwww = ZHwwHwww, 

HywHwwww = 2H V wwHwww, 
see |13| . These equations are in involution and can be solved in closed form. Up to the action of 
a natural equivalence group, this provides a complete list of integrable potentials. 

Theorem 1 [16] The generic integrable •potential of type (0|) is given by the formula 

H = Z(V + W) + eZ(V + eW) + e 2 Z(V + e 2 W) 

where e = e 2m ' 3 and Z"(s) = C( s )- He-re C, is the Weierstrass zeta- function, £' = — p, (p') 2 = 
4p 3 — g%. Degenerations of this solution correspond to 

H = {V + W) log(V + W) + e(V + eW) log(V + eW) + e 2 (V + e 2 W) log(V + e 2 W), 

H = ±V\(W), H=^, H = (V + W)ln(V + W), 
as well as the following polynomial potentials: 



ct, r * TT ,„„ 1. 



■', 



H = V Z W Z , H = VW Z + -W°, H = VW + -W 
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We refer to Sect. 2 for further details. 

In Sect. 3 we demonstrate that for systems of type ([lUp the integrability conditions take the 
form 

Hvvvv 



H 



vvvw 



H V v 

2Hyy\yHyyy 



H 



vv 



tt _ ZFfyyw f-, \ 

Hyy WW — — — , {12) 

Hyy 

2>Hy WW Hyy W - H WWW Hy V y 

HyWWW = 77 , 

Hyy 

tj ^>H VWW - 4H WWW Hy VW 
HWWWW = 



H 



vv 



These equations are also in involution, and can be solved in closed form. Up to the action of a 
natural equivalence group, we obtain a complete list of integrable potentials. 

Theorem 2 The generic integrable potential of type Mb]) is given by the formula 

V 



H = Vln 



a(W) 



where a is the Weierstrass sigma-function: a' jo = Q, Q' = — p, p' 2 = 4p 3 — (73. Its degenerations 
correspond to 

H = Vln^, H = VlnV, H = ^ + aW\ 

as well as the following polynomial potentials: 

V 2 VW 2 W 4 TT V 2 W 3 

H= 1 1 , H= 1 . 

2 2 4' 26 

Further details, as well as the dispersionless Lax pairs corresponding to integrable potentials 
from Theorem 2, are provided in Sect. 3. 

Remark. For H = V In V the equations ([TO]) take the form Vt = VW X , V x = —W v , implying the 
Boyer-Finley equations for V: V xx + (lnV)t y = 0. Similarly, the choice H = ^- +t- results in the 
equations V t = W x , V x = WW t — W y , implying the dKP equation for W: W xx = (WW t — W y )t- 
These Hamiltonian representations have appeared previously in the literature, see e.g. [3]. 

The case (jlip turns out to be considerably more complicated. The integrability conditions 
constitute an involutive system of fourth order PDEs for the potential H which is not presented 
here due to its complexity (see Sect. 4). However, it possesses a remarkable SL(3, i?)-invariance 
which reflects the invariance of the Hamiltonian formalism ([8]) under linear transformations of 
the independent variables x,y,t. Based on this invariance, we were able to classify integrable 
potentials. 

Theorem 3 The generic integrable potential H(V, W) of type \11\) is given by the series 

1 / 1 1 

" WglV) V 1 + gWW 6 + 92{V)W' 2 + " 

Here g(V) satisfies the fourth order ODE, 

g""(2gg 12 - g 2 g") + 2g 2 g'" 2 - 20gg'g" g >" + 16</V" + 18 99 " 3 - I8g' 2 g" 2 = 0, 
whose general solution can be represented in parametric form as 

VJ 2 (t) 

where w\(t) and 102 (t) are two linearly independent solutions of the hypergeometric equation 

2 

t(l - t)d 2 w/dt 2 - -w = 0. 
9 

The coefficients ffi(V) are certain explicit expressions in terms of g(V), e.g., g\ (V) = gg" — 2g' 2 
and so on. Degenerations of this solution correspond to 

H =Wg\v-y H = WV> H = V ~ W ^ W ' H = V-W 2 /2. 



We refer to Sect. 4 for further details. 

In Sect. 5 we provide a parametrisation of the generic integrable potential H(V, W) by 
generalised hypergeometric functions: 

Theorem 4 The generic integrable potential H(V, W) of type \11\) can be parametrised by gen- 
eralised hypergeometric functions, 

H = G(^\ V = £, W=£. (13) 

V u 2 - 1 / no ho 

Here G'(t) = 2 y 3 2 y 3 and ho,hi,h2 are three linearly independent solution of the hypergeo- 
metric system 

h U 4 h 2 h 2 «i(m-l) fc , 2 u 2 {u 2 -l) u 

n U 4 fc ^ 2«i(«l-l), , 2U2(«2-1), 

u 2 (l - U2)n U2 , a2 - -u 2 n U2 - -h = h Ul + ft U2 , 

3 9 3 u\ — ii2 3 U2 — ui 

^ Z tl U2 ""Hi 
3 U2 — U\ 
which can be viewed as a natural two-dimensional generalisation of the hypergeometric equation 
t(l-t)d 2 w/dt 2 - |w = 0. 

This parametrisation is based on a novel construction of integrable quasilinear systems in Go- 
dunov's form in terms of generalised hypergeometric functions which builds on [27, 29J. We 
believe that this construction is of independent interest, and can be applied to a whole variety 
of similar classification problems. The details are provided in Sect. 5 (Theorem 5). 

Our results lead to the following general observations: 
(i) The moduli spaces of integrable Hamiltonians in 2D are finite dimensional, 
(ii) In the two-component case, the actions of the natural equivalence groups on the moduli 
spaces of integrable Hamiltonians possess open orbits. This leads to a remarkable conclusion 
that for any two-component Poisson bracket in 2D there exists a unique 'generic' integrable 
Hamiltonian, all other cases can be obtained as its appropriate degenerations. 

We anticipate that the results of this paper will find applications in the theory of infinite- 
dimensional Frobenius manifolds, see j8] for the first steps in this direction. 



2 Hamiltonian systems of type ([6]) 

In this section we review the classification of integrable systems of the form 

/ v \ ( d/dx \ ( K 
\w) t ~\ d/dy ) V h w 

or, explicitly, 

vt = (hv)x, w t = (h w ) y . 



(14) 



The integrability conditions were first derived in |13| based on the method of hydrodynamic 
reductions. These conditions constitute a system of fourth order PDEs for the Hamiltonian 
density h(v,w), 

~\~otl vv il vw U vvw ^H vv ll ww il vvv il vvw <lvw'l"ww"< vvv i 
" J ww\" J wu} f^vv flivw ) ^vvvw — ^vw f^vvv {^vv f^www ~r i^ww^vvw ) 

-\-oH vw ti vvw ^ii vv n ww ii vvv ii vww -\- ti vw rtvvv 'T'vww ) 

h (h 2 —h h )h — Ah 2 h h 

ii>vw \ lt vw n vv ' b ww ) ">vvww — ^"'vw n vvw ' <"vww 

i^vv^ww \" / vwft"WWW ~r *^ww" j vww) '^ww'^vvv\'^vw'^vww t 'Ivv'lwwwji 

'T'vwy'lvw f^vv "'ww ) f^v www — "'vw'l'wwwy'l'ww'l'vvv ~T~ '^vv'^vww) 
M'ih 2 h 2 —Ihhh h A- h 2 h h 

<°"'VWVWW £>"'VV n 'WW"'WWW""OVW \ II ' VW "'WWW"'VVWJ 

'I'vwx'l'vw ""VV M ww ) 'T'ww ww — Q'lvw'T'Wwwy'l'Vw'I'Vww ^ww^vvw ) 

_4_^*/j h h 2 — ?h h h h — h h h 2 

^oiu ww iL vw ih vww i'%i)«(irai"tmi)tni%tmi) ll 'vw"'vv"'www 

This system is in involution, and its solution space is 10-dimensional. Eqs. (fl~4"|) can be repre- 
sented in compact form as 

sd h = d hds + 3h vw (2dh v dh w — h vw d h) det(dn) 

where s = h 2 w (h vv h ww — h 2 w ), d k h denotes A;-th symmetric differential of h, and n is the Hessian 
matrix of h. The contact symmetry group of the system (|14p is also 10-dimensional, consisting 
of the 8-parameter group of Lie-point symmetries, 

v — > av + b, 

w — >■ cw + d, 

h — ^ ah + j3v + ^w + 5, 

along with the two purely contact infinitesimal generators, 

d d d d d 

~^h v -~ &h w -~ h (fi — h v Q hv — h w Qh w )-Kj- + (fit; + h v Vl h )— — h (£l w + h w Slh)— — , 

with the generating functions Q = h 2 and ft = h 2 ^, respectively. It was observed in |13j that the 
integrability conditions simplify under the Legendre transformation, 

V = h v , W = h w , H = vh v + wh w — h, Hy = v, Hw = w, 

which brings Eqs. ([6]) into the Godunov form Q, 

(H v ) t = V x , (H w ) t = W y . 



The integrability conditions (|14p simplify to 

HvwHvvvv = ZHvvvHvvw, 

HvwHvvvw = ZHvvvHvww, 

HvwHvvww = HwwHyww + HvvvHwww, (15) 

HywHvwww = ZHvvwHwww, 

HywHwwww = IHywwHwww- 
This system is also in involution, and can be represented in compact form as 

Sd 4 H = d 3 HdS + 6H vw dVdW det{dN) 

where S = Hy W , and N is the Hessian matrix of H. The system (|15p is invariant under a 
10-parameter group of Lie-point symmetries, 

V -»• aV + b, 

W ^cW + d, 

H -»■ aH + pV 2 + ~fW 2 + fiV + vW + 5. 

Thus, both contact symmetry generators of the system (|14p are mapped by the Legendre trans- 
formation to point symmetries of the system (I15p . One can show that the action of the symmetry 
group on the moduli space of solutions of the system (|15p possesses an open orbit. The classifi- 
cation of integrable potentials H(V, W) is performed modulo this equivalence. 

2.1 Classification of integrable potentials 

The paper |16| provides a complete list of integrable potentials: 

Theorem 1 The generic integrable potential of type (OJ) is given by the formula 

H = Z{V + W) + eZ{V + eW) + e 2 Z(V + e 2 W) 

where e = e 2?n ' 3 and Z"(s) = C( s )- Here C, is the Weierstrass zeta-function, £' = — p, (p 1 ) 2 = 
4p 3 —53. Degenerations of this solution correspond to 

H = (V + W) log(V + W) + e(V + eW) log(F + eW) + e 2 {V + e 2 W) log(F + e 2 W), 

H=±V 2 ((W), H = ^, H = (V + W)\n(V + W), 
as well as the following polynomial potentials: 

H = V 2 W 2 , H = VW 2 + -W 5 , H = VW + -W 3 . 

5 6 

The corresponding systems Q possess dispersionless Lax pairs, see |16| for further details. 



3 Hamiltonian systems of type (|7J) 

In this section we consider Hamiltonian systems of the form (J2J) , 

/ v \ _ ( d/cfe \(K 

\ w ) t \ d/dx d/dy J \ h w 

or, explicitly, 

vt = (h w ) x , w t = (h v ) x + (h w ) y . 

Let us first mention two well-known examples which fall into this class. 
Example 1. The Hamiltonian 

w 2 



generates the system 

v t = w x , w t = e v v x + w y . 

Under the substitution v = u x , w = Ut these equations reduce to 



Utt ~ U t y 



e u q 



which is an equivalent form of the Boyer-Finley equation [6] 
Example 2. The Hamiltonian 

2^/2 



h = 1 wJw 

2 3 

generates the system 

1 1 

vt = —^=w x , w t = v x + 



^2w ' V2w y 

Introducing the variables V = V, W = y2w, we obtain 

v t = W x , WW t = V X + Wy. 

Setting W = Ut, V = u x we arrive at the dKP equation, 

u ty - u t u tt + u xx = 0. 

For the general system (J7|), the integrability conditions constitute a system of fourth order 
PDEs for the Hamiltonian density h(v,w), 

n>ww \>lvv ""WW 'T'vw)"'vvvv — Vtl vw ri vvw ' ^^ww^vv^vz 



WW">VV ,tj VVW 

_Od h h h _L h'- 

WW ""WWW K'VVV 



-\-Ah 2 h h —Ihhh h 4- h' z h 

~ L ±l b„„ n l l/wwil}' b WV *- jli WW n VW IL VVW n VVV I 'Vinin'" 



ft"ww y^vv^ww ^vw ) i^vvvw — *->fi vw ti vww ii vvw -\~ ori ww fi vww fi vvw ri vv 

3llww'lvw" j vvw ' " J vw^vvv'^www ~r ^ww f^vw ^v ww ^vvv ~r ^ww "'vvw ">vvv 5 
R'WW\P'vv""W'U) ""ww ) thrown) — ^'^•vw'^vww ~r *'"'vv""vww'^ww 

~Oll ww ri vw ll vww ll vvw -\~ ll ww ll vv ll www ll wvv -\~ tiwW'^VW'^WWW'^VVV T ^yjyj" J VWW" J VVV-! 

h (h h — h 2 )h — —1h 2 h h — 1h h 2 h 

ll'WW\ ll 'VV l <"WW ' b VW I "•VWVJW ^ n "UW IL WWW n 'VWW ' JII 'WW II 'VWW"'VW 

~rO/l ww il vv il www iT vww T" "ww" j vw'^ j www'^vvw t ^ww ""www l^vvv j 
h (h h — h 2 )h — -2h 2 h 2 — 2h h h h 

•l'WW\ , ''VV ll 'WW lh VW I ""WWWW Z '"'Vw"'WWW z '"'WW"'WWW ll 'VW"'VWW 

— 'Ih 2 h 2 -4- 3ft ft h 2 4- 4ft 2 ft ft 

°"'Ww"'VWW ' K>II 'WW"'VV""WWW ' ^ n ww "" www lb vv w 

10 



(16) 



As an illustration of the method of hydrodynamic reductions, details of the derivation of these 
conditions are presented in the Appendix. We have verified that system (|16p is in involution, and 
its solution space is 10-dimensional. The integrability conditions can be represented in compact 
form as 

sd h = d hds + 3h ww (2(dh w ) — h ww d h) det(dn) 

where s = h% mo {h vv hww — h% w ), d h denotes fe-th symmetric differential of h, and n is the Hessian 
matrix of h. This system is invariant under an 11-dimensional group of contact symmetries which 
consists of the 9-dimensional group of Lie-point symmetries, 

v — > av + b, 

w — )• pv + cw + d, 

h — ^ ah + j3v + 7U> + 5, 

along with the two purely contact infinitesimal generators, 

d d d d d 

~^hva ^h w ^- + (H - h v Q hv - h w Q hw )— + (Cl v + h v Q h )-^— + (Q w + h w VL h ) 



•dv nw dw v " nv w nw 'dh y v " nj dh v K w w n 'dhJ 

with the generating functions Vl = h v h w and Q = h"^, respectively. To solve this system we apply 
the Legendre transformation, 

V = h v , W = h w , H = vh v + wh w — h, Hy = v, H\y = w. 

The transformed equations ([7|) take the Godunov form (|10p . 

(H v ) t = W x , (H w ) t = V x + W y , 

and the integrability conditions (I16p simplify to 



Hvvw '- 
Hwvw 



^Hyyy 
Hyy 
2HyywHyyy 



H< 



vv 



tt _ zn VVW i-\ 7 \ 

tiyy 

3Hy WW Hyy W - H WWW Hyyy 

UVWWW = 7} j 

tiyy 

rr ® H VWW - 4H WWW Hyy W 

UWWWW = 77 • 

tiyy 

This system is also in involution, and can be represented in compact form as 

Sd A H = d 3 HdS - 6H vv (dW) 2 det(dN), 

where S = Hyy, and N is the Hessian matrix of H. The system (|17p is invariant under an 
11-parameter group of Lie-point symmetries, 

V -^aV + bW + p, 

W^cW + d, (18) 

H ^aH + (3W 2 + -fVW + eV + kW + rj. 

Thus, both contact symmetry generators of the system (|16p are mapped to point symmetries of 
the system (I17p . All our classification results will be performed modulo this equivalence. 

11 



3.1 Classification of integrable potentials 

The main result of this section is a complete classification of solutions of the system (|17p : 

Theorem 2 Modulo the equivalence group U8\) , the generic integrable potential H(V, W) is given 
by the formula 

where a is the Weierstrass sigma-function: a' jo = £, £' = — p, p' 2 = 4p 3 — 53. Its degenerations 
correspond to 

H = Vln^, (20) 

H = VlnV, (21) 

H=^ + aW*, (22) 

as we// as the following polynomial potentials: 

v 2 vw 2 w 4 

H = — + —— + — , (23) 

2 2 4 v ' 

V 2 W 3 
*= T + -jj- (24) 

Proof: 

In the analysis of Eqs. (|17p it is convenient to consider two cases, Hyyy 7^ and Hyyy = 0. 
Case 1: Hyyy 7^ 0. Then Eqs. (fTTjh and (fTT)) ? imply -ffyyv' = cH vv , c = const, so that 
ifyy = — M7T77WT- The substitution into (|17p ^i implies that /(W) must be linear. Modulo the 
equivalence group we can thus assume that Hyy = y, which gives 

H = V \n(V) + g(W)V + h(W). 

The substitution of this ansatz into {\17$ a gives h'" = 0, so that we can set h equal to zero modulo 
the equivalence group. Ultimately, the substitution into the last equation (|17p^ gives 

//// n "2 

9 =% . 

The general solution of this equation is g = —lna(W), where a is the Weierstrass sigma- 
function defined as (lno - )' = £, £' = — p with p satisfying p' 2 = 4p 3 — 53 (notice that 52 = 0). 
Degenerations of this solution give g = — In W and <7 = 0. This leads to the three cases 

H = Vln-^—, H = Vln^-, H = VlnV. 
Case 2: Hyyy = 0. Then Hyy = a(W), and (fl7> implies 

a = 2 — , 



12 



so that a{W) = w+ . There are again two subcases to consider depending on whether Hyvw 
is zero or non-zero. If Hyvw 7^ then p^ and, modulo the equivalence group, one can set 
Hyy = ttf. Integrating twice with respect to V we obtain 

H=^ + b(W)V + c{W). 

The substitution into <\17\ a gives b'" = —3b"/W. Modulo the equivalence group, one can set 

b(W) = l/W. Thus we get 

V 2 V 
H = h — + c(W). 

By a translation of V one can remove the intermediate term X, leaving H = Xtf + c(W). Then 
the substitution into (|17p ^ gives 

Ac'" 

w 
Modulo the equivalence group, this gives 

H = w + aW7 - 

In the second subcase, Hyyy = Hyyyy = 0, one can normalize H as H = V 2 /2 + a{W)V + b(W) . 
With this ansatz Eqs. (|17|) i - (|17p ^ are trivial, while the last two equations imply, modulo the 
equivalence group, that 

V 2 
H= — + aVW 2 + a 2 W 4 + (3W 3 . 



The case a = leads to 



V 2 W 3 

H = 1 , 

2 6 



The case a^O leads, on appropriate rescalings, to 

TT V 2 VW 2 W 4 
2 2 4 

This finishes the proof of Theorem 2. 

3.2 Dispersionless Lax pairs 

We recall that a quasilinear system ([TJ is said to possess a dispersionless Lax pair, 

S x = F{u,S t ), S y = G(u,S t ), (25) 

if it can be recovered from the consistency condition S xy = S yx (we point out that the depen- 
dence of F and G on St is generally non-linear). Dispersionless Lax pairs first appeared in the 
construction of the universal Whitham hierarchy, see |23| and references therein. It was observed 
in |39j that Lax pairs of this kind naturally arise from the usual 'solitonic' Lax pairs in the 
dispersionless limit. It was demonstrated in |14| [15] that, for a number of particularly interest- 
ing classes of systems, the existence of a dispersionless Lax pair is equivalent to the existence of 
hydrodynamic reductions and, thus, to the integrability. Lax pairs form a basis of the dispersion- 
less 9-dressing method [5j and a novel version of the inverse scattering transform associated with 
pararmeter-dependent vector fields |24j. In this section we calculate dispersionless Lax pairs, 

S x = F{V, W, S t ), S y = G(V, W, S t ), 
13 



for all integrable potentials appearing in Theorem 2. Thus, we require that the consistency 
condition S xy = S yx results in the corresponding equations (|10p . Without going into details of 
calculations, which are quite standard, we only state final results. Let us point out that in all 
examples discussed below the first equation of the Lax pair does not explicitly depend on V. 
The general class of Lax pairs of this kind was discussed in |28| . 

Case flSH): H = ^- + ^-. The equations (JTDJ take the form 

V t = W x , 

WW t = V X +Wy. 

The corresponding dispersionless Lax pair is 



o2 

s x = f-w, 

Sy = ~^+S t W + V. 



Case 0235: H = \ + ^4f- + \. The equations ([ID]) take the form 

Vt + WW t = w x 

VW t + WV t + 3W 2 W t = V x + W y . 

The corresponding dispersionless Lax pair is 

s x = -f-s t w, 



o7 

S y = -^r + WSt + St(V-W 2 ) 



Case (1221) : H = ^, a = 0. The equations (fTOjl take the form 



V t VW t 



w„ 



X: 



w 2 w 3 

The corresponding dispersionless Lax pair is 



W W 2 
vv t V 2 W t rr „ r 



w 2 

Sx ~ ~2S?' 

vw w 5 

y ~ ~w + w 

Case (1221) : H = ^w + aW 7 . Without any loss of generality we will set a = 1/168. The equations 



2W 
TUl) take the form 



v t vw t 



w w 2 



w x , 



-5 



Wt V 2 Wt 1 
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The corresponding dispersionless Lax pair is 

aW 2 



b x 



2 ' 



aa'W 5 



S v = aVW m . 

here a = a(St) satisfies the ODE 3a' 2 — 5 = 2aa". 

Case d2T]): H = VlogV". The equations JTD]) take the form 

v t = vw x , 

= V X + Wy. 

The corresponding dispersionless Lax pair is 

S x = -\og{W + S t ), 
S - _J1_ 

v ~ w + s t - 



Case fl20|): H = V log ^. The equations (JTDJ) take the form 






it 



.T- 



-vd~4 = v * +w *- 

The corresponding dispersionless Lax pair is 

S x = -ln(S t + W)-eln(S t + eW)-e\n(S t + e 2 W), e = e^ , 

_ vws t 

Case ()19|) : H = Flog a Xyy We recall that (log a)' = ( and (' = — p where p = p(w,0,g^) is 
the Weierstrass p- function, p' 2 = 4p 3 — g%. The equations fjlUp take the form 

^V t -((W)W t = w x , 
-C(W)V t + Vp(W)W t = v x + w y . 

The corresponding dispersionless Lax pair is of the form 

S x = F(W,St), 
S y = VA(W,St), 

where the functions F and A satisfy the equations 

F W = -A, Fs = ^-((W) 



and 



^ = pW-f, A ww = 2^f-2Ap(W), 



15 



respectively (here £ = St). Setting A = 1/u one can rewrite the last two equations for A in the 
equivalent form, 

«£ = uw - p(W)u 2 , u W w = 2p(W)u. (26) 

Here the second equation is a particular case of the Lame equation. It has two linearly indepen- 
dent solutions |3], 



a{W)a(ay a{W)a{a)' 

where a is the zero of the p-function: p(a) = 0. Thus, one can set 



<r(V^)ff(a) v ' <r(W>(a) 

The substitution into (|26j) i implies a pair of ODEs for a(£) and b(£), 

a' = p'(a)6 2 , 6' = -p'(a)a 2 . 
This system can be solved in the form 



6(0- 



Tp'(£,0,g 3 )-1 



7p(£,0,# 3 ) 






>P'(£,0,53) + 
where the constants c, r and 7 are defined as 



T = Ts 9 ~ s 



1/2 



7 = 25. 



Setting 53 = 1, p'{a) = i one obtains 



«(0 



,^ p'(g,0,l)->/3 



P 



(£,0,1) + ^' 



-1/3 
3 ' 



KO 



V^ 1/6 

p'(a)" 3 ' 



P 



6p(g,0,l) 
(£,0,l) + v / 3' 



Ultimately, 



u = 1 



/g c -WC(a) ^(W - + a) g/(fl - \/3 _ 6 . c h/c(«) g (W ~ g) P(fl 

a(W>(a) p'(£) + v^ ff(W>(a) g/(0 + \/3" 



The functions .A and i 7 can be reconstructed via A = 1/u, 



u w 



F w = -l/u, Ft = -^-C(W). 
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4 Hamiltonian systems of type ([8]) 

In this section we consider Hamiltonian systems of the form (JSJ) 



v 
w 



2b wW / « \ _^_ 

u> J dx \ v 2w J dy \ w x w y 



or, explicitly, 



h v 

il,n 



v t = (2vh v + why, - h) x + (vh w ) y , w t = (wh v ) x + (2wh w + u/i„ - h) y . 
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The integrability conditions constitute a system of fourth order PDEs for the Hamiltonian density 
h(v, w) which is not presented here due to its complexity. We have verified that this system is in 
involution, and its solution space is 10- dimensional. It is invariant under an 8-dimensional group 
of Lie-point symmetries, 

v — > av + bw, 
w — >■ cv + dw, 

h — ^ ah + j3v + ^w + 5, 

along with the two purely contact infinitesimal generators, 

d d d d d 

-0, hv - ^h w -z— + (ft - h v Q hv - h w 9, hw )-— + (ft„ + h v tt h )-^— + (Q w + h w tt h ) 



•dv riw dw ' y v Hv w nwJ dh ' v v ' "° n 'dK y w w n 'dhJ 

with the generating functions ft = h v (vh v + wh w ) and ft = h w (vh v + wh w ), respectively. As in 
the previous two cases, the integrability conditions for the Hamiltonian density h(v, w) simplify 
under the Legendre transformation, 

V = h v , W = h w , H = vh v + wh w — h, Hy = v, H\y = w, 

which brings the corresponding equations (JBJ) into the Godunov form (HID . 

(H v ) t = (VH V + H) x + (WH v ) y , (H w ) t = (VH W ) X + (WH W + H) y . 

Remark. Equations (|lip coincide with the so-called EPDiff equations [22j, 

rri( — u x curl m + V(u, m) + m divu = 0, 

where u = — (V, W, 0) and m = (Hy,H\y,0). These equations have a geometric interpretation 
as the Euler-Poincare equations for the geodesic motion on the group of diffeomorphisms. 
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The integrability conditions for the system (|11|) take the form 
SHwvv =2(HwHwv — HyHyyw) + 4H V (H VVW — HyyyHyww) 

+ Y2H\\/HyyHyyy/ — 12H\y Hyy Hyw Hyyy — 2AHy Hyy Hy\y Hyy\y 
+ 8HyHy W Hyyy + 6HyHyyHywW + WHy Hyy H]y\y Hyyy 

+ 12Hyy(Hy W - H VV H WW ), 

SHyyyw =Hyyy(2H w Hyyw ~ AHyHy/HywW ~ HyHy/Ww) + ^HyHyyw HywW 

— ZHwHyyHyy/Hyyw ~ ^H\yHy W Hyyy — 6HyHy W Hyyw 

15 3 

+ -^-H\yHyyHywW — QHyHyyHywHywW + ^HyHyyHwWW 

1 3 

~ 7} H\V Hyy H\yw Hyyy + -HyHyyHwwHyyw + 9Hy Hyw Hww Hyyy 

+ 12HyyHy W (H vw - HyyH W w), 

SHyywW =2{HyHywW — H\yHyy\y) + 2HyH\y{Hyy\yHy\y\y — HyyyHwWw) 

— 8Hy W (H\yHyy\y + HyHyww) + 6H\y Hyy Hyw HwWV + QHyH\y\yHy\yHyy\y 

— HyyHww(H\yHyy\y + HyHyww) + 3HwHyyH\yWW + 3HyH ww Hyyy 
+ ^{2Hy W + Hyy Hww)(Hy W — HyyHww), 

SHyy/WW =HwWw{2HyHyy/W ~ 4:HyH\yHyy\y ~ H w Hyyy) + 3H w Hyy\yHy\yW 

— ZHyHwwHywHywW ~ ^HyHy W Hyi/WW ~ 6H\yHy W Hyww 
15 3 

+ -yHyH ww Hyy W - 6H W H WW Hy W Hyy W + -H W H ww Hyyy 

1 3 

— -HyHyyHwwHwWW + ^H\y Hyy Hy/W HvWW + 9HwHy\yHyyHwWW 
+ 12HwwHyw{H vw — HyyHww), 

SH WWWW =2(H V H WWW - H W H VWW ) + 4H W (H VWW - H www Hy VW ) 

+ 12HyH ww Hy\yw — 12HyH\ywHywHwWW — 24H\yH\ywHywHywW 
+ 8H\yHy W H\yWW + 6H\yH ww Hyy\y + 10H\yH\ywHyyH\yWW 



+ 12H WW (H VW — HyyHww), 



where 



S = H w Hyy — 2HyHwHyw + HyHww 



This system is manifestly symmetric under the interchange of V and W, and can be represented 
in compact form as 



Sd A H = 2d 3 HdS - 6(dH) 2 det(diV) + 6dHd 2 Hd(det N)+ 
12(H v dH w - H w dHy)(d 2 H v dHw - d 2 H w dH v ) - 12det N(d 2 H) 2 , 

18 



(27) 



where N is the Hessian matrix of H. We have verified that the system (|27p is in involution, 
and is invariant under a 10-dimensional group of Lie-point symmetries generated by projective 
transformations of V and W along with affine transformations of H, 

aV + bW + c aV + (3W + 7 

V ->■ — — , W -> — — , H —■ [iH + v. 

pV + qW + r pV + qW + r 

The corresponding infinitesimal generators are: 

d d 

2 translations : 



d d d d 

4 linear transformations : V——, W——, V——, W- 



dV' 9V" dW 8W' 

d d d d 

2 projective transformations : V —— + VW——, VW-— + W 7^777; 

oV oW oV oW 



(28) 



2 affine transformations of H : 7—-, H——. 

on oH 

Projective transformations (|28|) constitute the group of 'canonical' transformations of equations 
()lip : combined with appropriate linear changes of the independent variables x,y,t, they leave 
equations (|lip form-invariant. Let us consider, for instance, the projective transformation V = 
1/V, W = W/V. A direct calculation shows that the transformed equations take the form 

(VHy + H) t = (Hy) x + (WHy) y , (VH w ) t = {H#) x + (W H w + H) y . 

They assume the original form (|lip on the identification t — >• x, x — >• t, y — )■ —y. In other words, 
the projective invariance of the integrability conditions is a manifestation of the invariance of the 
Hamiltonian formalism ([8]) under arbitrary linear transformations of the independent variables. 
This is analogous to the well-known invariance of Hamiltonian structures of hydrodynamic type 
in 1 + 1 dimensions under linear transformations of x and t |34] , We emphasise that this invariance 
is not present in the case of constant coefficient Poisson brackets. 

4.1 Classification of integrable potentials 

The analysis of this section is somewhat similar to the classification of integrable Lagrangians 
of the form J utg(u x ,u y ) dxdydt proposed in [17] , This suggests that there may be a closer link 
between these two classes of equations. The main result of this section is the following 

Theorem 3 The generic integrable potential of type \11\) is given by the series 

H = WgJV) V 1 + 9l (V)W Q + g 2 (V)W^ + " 
Here g(V) satisfies the fourth order ODE, 

g""(2gg' 2 - g 2 g") + 2g 2 g'" 2 - 20gg'g"g'" + 16g> 3 g>" + 18^" 3 - 18g' 2 g" 2 
whose general solution can be represented in parametric form as 

w 2 {t) 
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where wi(t) and W2(t) are two linearly independent solutions of the hypergeometric equation 
t(l — t)d 2 w/dt 2 — gW = 0. The coefficients giiV) are certain explicit expressions in terms of 
giV), e.g., giiV) = gg" — 2g' 2 and so on. Degenerations of this solution correspond to 

H =wgly-y H = Wv> H = v~wio g w, H = V-W 2 /2. 

In the rest of this section we provide details of the classification, and discuss various properties 
and representations of the generic solution. The system (|27p for H(V, W) is not straightforward 
to solve explicitly. We will start with the investigation of special solutions which are invariant 
under various one-parameter subgroups of the equivalence group. In this case the system of 
integrability conditions for H(V, W) reduces to ODEs which are easier to solve. Up to conjugation 
and normalisation, there exist four essentially different one-parameter subgroups of the projective 
group SL(3), with the infinitesimal generators 

Tr d TTT d _ d d TTT d d d 

aV h W V 1 W 1 

dV dW : dV dW 8V dW dV 

Combined with the operators d/dH, Hd/dH this leads to the following list of eleven essentially 
different 'ansatzes' governing invariant solutions (in what follows we do not consider degenerate 
solutions for which the expression S = H^Hyy — 2HyH]yHy]y + ^%w equals zero): 
Case 1. Solutions invariant under the operator d/dV + d/dH are described by the ansatz 
H = V + F(W). The integrability conditions imply F"F"" — 2F'" 2 = 0. Modulo the equivalence 
transformations this leads to integrable potentials H = V — W 2 /2 and H = V — W log W, which 
constitute the last two cases of Theorem 3. Applying to the first potential transformations from 
the equivalence group we obtain integrable potentials of the form 

H(V W) = Q ^ W) 

{ ' J l 2 (V,W) 

where Q and I are quadratic and linear forms, respectively (not necessarily homogeneous). The 

integrability implies that the line I = is tangential to the conic Q = on the V, W plane. Any 

such potential can be reduced to the form H = V — W 2 /2 by a projective transformation which 

sends I to the line at infinity. 

Case 2. Solutions invariant under the operator d/dW + Hd/dH are described by the ansatz 

H = e w FiV). This case gives no non-trivial solutions. 

Case 3. Solutions invariant under the operator Wd/dV + d/dW are described by the ansatz 

H = F(V — W 2 /2). A simple analysis leads to the only polynomial potential H = V — W 2 /2, 

the same as in Case 1. 

Case 4. Solutions invariant under the operator Wd/dV + d/dW + d/dH are described by the 

ansatz H = W + F(V — W 2 /2). In this case F turns out to be linear so that, modulo translations 

in W, we again arrive at the same potential as in Case 1. 

Case 5. Solutions invariant under the operator Wd/dV + d/dW + Hd/dH are described by 

the ansatz H = e F(V — W 2 /2). A detailed analysis shows that this case gives no non-trivial 

solutions. 

Case 6. Solutions invariant under the operator Vd/dV + d/dW are described by the ansatz 

H = F(W — In V). This case gives no non-trivial solutions. 

Case 7. Solutions invariant under the operator Vd/dV + d/dW + d/dH are described by the 

ansatz H = W + F(W — lnU). This case gives no non-trivial solutions. 
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Case 8. Solutions invariant under the operator Vd/dV + d/dW + [iHd/dH are described by 

the ansatz H = e^ F(W — InV), also no non-trivial solutions. 

Case 9. Solutions invariant under the operator aVd/dV + Wd/dW are described by the 

ansatz H = F(W a /V). Here one can assume a/0,1. A straightforward substitution implies 

that, without any loss of generality, one can assume F to be linear, while the parameter a can 

only take two values: a = 2 or a = — 1. This results in the two rational potentials H = W 2 /V 

and H = 1/VW. Notice that they are related by the projective transformation W —> 1/W, V — > 

V/W. 

Case 10. Solutions invariant under the operator aVd/dV + Wd/dW + d/dH are described 

by the ansatz H = \n.W + F(W a /V). This case gives no non-trivial solutions. 

Case 11. Solutions invariant under the operator aVd/dV + Wd/dW + /j,Hd/dH are described 

by the ansatz H = W^F{W a /V). A detailed analysis shows that, modulo the equivalence group 

and the solutions already discussed above, the only essentially new possibility corresponds to the 

ansatz H = w ,y, (here /i = —1, a = 0). It leads to the fourth order ODE for g = g(z), 

g""(2gg' 2 - g 2 g") + 2g 2 g'" 2 - 20gg'g"g'" + 16g' 3 g'" + 18gg" 3 - 18g' 2 g" 2 = 0, (29) 

which possesses a remarkable SL(2, i?)-invariance inherited from 

az + /3 



"fz + 5 



~g = (~iz + S)g; (30) 



here a, /3,7, 5 are arbitrary constants such that a5 — (3j = 1. Moreover, there is an obvious 
scaling symmetry g — > Xg. The equation (|29p can be linearised as follows. Introducing h = g' /g, 
which means factoring out the scaling symmetry, we first rewrite it in the form 

ti"(ti - h 2 ) - 2h" 2 + Uhh'h" - Ah 3 h" - 15/i' 3 + 9h 2 ti 2 - 3h 4 ti + h e = 0; (31) 

the corresponding symmetry group modifies to 

az + j3 



721 + 5 



h={ 1 z + 5Yh + 1 { 1 z + 8). (32) 



We point out that the same symmetry occurs in the case of the Chazy equation, see [T], p. 342. 
The presence of the S'L(2,i?)-symmetry of this type implies the linearisability of the equation 
under study. One can formulate the following general statement which is, in fact, contained in 

m. 

Proposition 1. Any third order ODE of the form F(z,h,h' ,h" ,h'") = 0, which is invariant 
under the action of SL(2, R) as specified by &32\) . can be linearised by a substitution 

h = —lnu)2, z = — (33) 

dz tt>2 

where u>i(t) and W2(t) are two linearly independent solutions of a linear equation d 2 w/dt 2 = 
V(t)w with the Wronskian W normalised as W = W2dw\ / dt — W\dw2 / dt = 1. The potential V{t) 
depends on the given third order ODE, and can be efficiently reconstructed. 

In particular, the general solution of the equation \31\) is given by parametric formulae h33\) 
where u>i(i) and ui2{t) are two linearly independent solutions of the hypergeometric equation 
d 2 w/dt 2 = ^ t ri_ t \ W with W = 1. 
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Proof: 

Our presentation follows |17| . Let us consider a linear ODE d 2 w/dt 2 = V(t)w, take two linearly 
independent solutions wi(t), W2(t) with the Wronskian W = 1, and introduce new dependent 
and independent variables h, z by parametric relations 

u d i Wl 

h = — — In W2 , z = — . 

dz ' W2 

Using the formulae dt/dz = w 2 and h = W2dw2/dt, one obtains the identities 

ti -h 2 = w\ V, 

h" - 6hti + Ah 3 = wl dV/dt, 
h!" - I2hh" - G(ti) 2 + 48h 2 ti - 24/i 4 = w% d 2 V/dt 2 , 
where prime denotes differentiation with respect to z. Thus, one arrives at the relations 

_ (h" - 6hti + Ah 3 ) 2 _ (dV/dt) 2 
1 ~ [h! - h 2 f ~ V 3 ' 

_ h'" - Ylhh" - 6(ti) 2 + A8h 2 h' - 24/i 4 _ d 2 V/dt 2 
2 ~ [h! - h 2 ) 2 ~ V 2 ' 

We point out that I\ and I2 are the simplest second and third order differential invariants of the 
action (|32p whose infinitesimal generators, prolonged to the third order jets z,h,h' ,h",h"' , are 
of the form 

Xi=d g , X 2 = zd z - hd h - 2tid h > - 3h"d h » - Ah'"d h »,, 

X 3 = z 2 d z - (2zh + l)d h - (2/j + Azh')dhi - (6ti + §zh")d h „ - (12h" + »zh m )d h >»\ 

notice the standard commutation relations [Xi,^] = Xi, [Xi,Xg] = 2X2, [.X^-X^] = X3. One 

can verify that the Lie derivatives of I\, I2 with respect to X\, X2,X^ are indeed zero. Thus, any 

third order ODE which is invariant under the S'L(2,i?)-action (|32p . can be represented in the 

form I2 = F(I\) where F is an arbitrary function of one variable. The corresponding potential 

t^m ■* 1 ■ d 2 V/dt 2 ^{(dV/dt) 2 * 
V (t) has to satisfy the equation k — = r ' 



v 2 V v 3 

This simple scheme produces many the well-known equations, for instance, the relation I2 = 

—24 implies the Chazy equation for h, that is, h'" — 12hh" + 18(/i') 2 = 0. The corresponding 

potential satisfies the equation d 2 V/dt 2 = —2AV 2 . 

(h" - 2hh') 2 

Similarly, the choice I 2 = h - 8 results in the ODE h'" = Ahh" - 2(ti) 2 + ^— ~^- 

h! — h 2 
which, under the substitution h = y/2, coincides with the equation (4.7) from [2]. The potential 

V satisfies the equation Vd 2 V/dt 2 = (dV/dt) 2 - 8V 3 . 

The relation I 2 = h -9 gives ti"(ti-h 2 ) = h 6 -3h 4 h' +9h 2 (ti) 2 -3(h'f -4h 3 h" + {ti') 2 . This 
equation appeared in |17| in the context of first order integrable Lagrangians. The corresponding 
potential V satisfies the equation Vd 2 V/dt 2 = (dV/dt) 2 - W 3 . 

Finally, the relation I2 = 2I\ + 9 coincides with the equation (J3T|) . The corresponding 
potential V satisfies the equation Vd 2 V/dt 2 = 2(dV/dt) 2 + 9V 3 . It remains to point out that 
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the general solution to the last equation for V is given by V = —9 .2 , t+b - Without any loss 

2 1 

9 t(l-t) ■ 



of generality one can set V = q f n_ f \ ■ It this case the linear equation d 2 w/dt 2 = V(t)w takes 



the hypergeometric form corresponding to the parameter values a = — tt, 6=—^, c = 0: 
t(l — t)d 2 w/dt 2 — gW = 0. This finishes the proof of Proposition 1. 

As h = g' /g, this immediately implies the following formula for the general solution of 
Proposition 2. The general solution of the equation A29\) is given by parametric formulae 

W\ 

g = w 2 , z = — , 
w 2 



where W\ and w 2 are two linearly independent solutions to the hypergeometric equation t(l 
t)d 2 w/dt 2 - § 10 = 0. 

Remark. Expansions at zero give 



9 t 2 211t 3 

= w<2 In t H 1 1 h • • • , wo 

2 3 1458 


t 2 10t 3 

" t+ 9 + 243 + -' 


t 2 10t 3 wi , 9 

1 h ..., = — = ln*H 

9 243 ' w 2 2t 


1 lit 34t 2 715t 3 

2 + 54 ' 729 ' 39366 



W\ 

so that 

q = W2 = t H h — h ..., 2; = — - = hit + - 1 h h — h .... 

y 9 243 ' w 2 2*254729 39366 

Solving the first relation for t in terms of g and substituting into the second, one gets an implicit 
relation connecting g and z, 

9 1 112o 289a 2 4381o 3 

z = Inq 1 1 1 1 h ■ • ■ • 

y g 2 27 486 39366 

Similarly, expansions at infinity give 

«*=Wi-i-4*--.. > ), w 2 =t^d-l 



3* 45t 2 "7 ' V 6t 126t 2 

Hl = t i/s( 1 _l 41 



?<>2 



6t 1260t 2 



so that g = w 2 can be represented explicitly as 



, a 1165 a 2 5280035 a 3 

g(z) = z 1 - — 



143 z 12 46189 z 18 

here a = 1/140. For other values of a this formula represents the general solution to (|29p in the 
form g = z(\ — a/z 6 — b/z 12 — c/z 18 — ...). The corresponding potential H(V, W) takes the form 

It can be viewed as a perturbation of the potential H = yyj constructed before. Computer 
experiments show that the 'generic' integrable potential H(V, W) can be represented as 

TT v^ a jk 1/11 24 1308 1308 



j,fc>0 


FW 


/ 1 1 

V f 6 + w 6 ' 
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An alternative representation of this solution, 

1 / 1 1 

" W^V) V 1 + 91 (V)WG + g 2 (V)W^ + "' 

can be obtained be rearranging terms in the above sum. The substitution of this expression into 
the integrability conditions implies that g(V) has to satisfy the fourth order ODE (|29|) . g\{V) 
is expressed in terms of g(V) via the Rankin-Cohen-type operation, gi(V) = gg" — 2g' 2 (recall 
that, due to (|30p . g(V) transforms as a modular form of weight 1), and so on. This finishes the 
proof of Theorem 3. 

In Sect. 5 we provide a parametrisation of the generic integrable potential by generalised 
hypergeometric functions. 

4.2 Dispersionless Lax pairs 

Here we present Lax pairs for some of the simplest potentials found in the previous section. 
Case H = V — W 2 /2. The corresponding system (fTTj) takes the form 

W y = WW X - 2V X , W t = VW X - 5WV X + 3W 2 W X - V y , 

and possesses the Lax pair 

S y = -WS x + ± S t = {V-W 2 )S x + ^-^- E . 

Case H = V — W\ogW. The corresponding system (|11|) takes the form 

W y = (WhiW-2V) x , W t /W = (VlnW + V) x + {2WlnW + W -V) y , 
and possesses the Lax pair 

S y = -S x (lnW + l)+p(S x ), S t = S X (V -W\nW-W)- ^S xP '(S x )W 

where p(z) satisfies the ODE zp" + p' 2 + 3p' = 0. This gives p' = a _. so that 

p(S x ) = ln(S x -l) + eln(S x -e)+e 2 ln(S x -e 2 ), e = eT. (34) 

Case H = r^y- The corresponding system (|11|) takes the form 



1 \ ( 1 \ / 1 



. 



v 2 wj t \v 2 ) y ' \vw 2 ) t \W 2 JJ 

and possesses the Lax pair 

S x = a(S t )/V, S y = b(S t )/W 

where the functions a{z) and b(z), z = St, satisfy a pair of ODEs 

a' = l-4, 6' = 1-A. 
26' 2a 

These equations can be solved in parametric form as 

a ( z ) - P'(P) + X h( z\ - P ' {P) - X z - -2CM 
Q{) 2p(p) ' b[Z) ~ 2p(p) ' Z ~ 2aPh 
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(35) 



where p(p) and C(p) are the Weierstrass functions, p' 2 = 4p 3 + A 2 , Q' = —p. The resulting Lax 
pair can be written in parametric form as 



or, equivalently, 



_ 1 p'jp) + A _ 1 p'{p) - A 



S x S y (VS x -WS y ) = ^, S x S y = ^, S t = -2((p). 



Notice that equations (j35j) coincide with the Euler-Lagrange equations corresponding to the 
Lagrangian density J % dxdydt upon setting V = #*-, W = #*- (we thank Maxim Pavlov for 

pointing this out). In this context, the above Lax pair appeared previously in |33j . 

5 Godunov systems and generalized hypergeometric functions 



In this section we develop the general theory of Godunov's systems [21J . and describe the Go- 
dunov form of n-component quasilinear systems constructed in |27[ [29] in terms of generalized 
hypergeometric functions. In particular, this provides a Godunov representation for any generic 
2-component integrable system |14[ 127] . Applied to the system (|lip . this construction gives a 
parametrisation of the generic integrable potential H(V, W) in the form fp~3|> . 

In Sect. 5.1 we recall the main aspects of the Godunov representation, and clarify its sym- 
metry properties. In Sect. 5.2 we construct the Godunov form for quasilinear systems found 
in |27|, 129] in terms of generalised hypergeometric functions. In Sect. 5.3 we specialise this 
construction to integrable potentials H(V, W) of the type (jlip . thus proving Theorem 4. 

5.1 Godunov systems 

A 2 + 1 dimensional quasilinear system in n unknowns v = (v 1 ,...,v n ) is said to possess a 
Godunov representation |21| if it can be written in the conservative form, 

(F 0ti )t + (F 14 ) x + (F 2ti )y = 0, i = l,...,n, (36) 

where potentials Fq, F\ , F 2 are functions of v, and F a j = dFa/dv 3 , a = 0, 1, 2. Any such system 
automatically possesses an extra conservation law, 

C(F ) t + £(Fi) x + C(F 2 ) y = 0, (37) 

where C denotes the Legendre transform: C(F a ) = F a }~v — F a . Many systems of physical 
origin are known to be representable in the Godunov form. This representation is widely used 
for analytical/numerical treatment of quasilinear systems. 

Note that the systems ([9|)- (|lip are written in the Godunov form. For example, in the case 
(fTTjl we have n = 2, v 1 = V, v 2 = W, F = -H, F 1 = VH, F 2 = WH. Recall that a 2 + 1 
dimensional quasilinear system of n equations for n unknowns possesses a Godunov form iff it 
possesses n+1 conservation laws of hydrodynamic type. The following fact will be useful: 

Proposition 1. The Godunov representation W6\) is form-invariant under the projective action 
of GL n+ \ defined as 



a 



Z(v)' ta l(yy 
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here l l , I are linear (inhomogeneous) forms in v. 

Proof: 

With any Godunov system we associate the following geometric objects. Let us consider an 
auxiliary (n + l)-dimensional affine space A n+1 with coordinates xi, . . . ,x n ,x n +i. With any 
potential F a we associate an n-parameter family of hyperplanes, 

x n+1 -x k v k + F a (v) = 0. (38) 

The envelope of this family is a hypersurface Mp C A n+1 defined parametrically as 

[X\, • • • , X ni X n ^i) — \-f(x,lj ■ • • ; ^a,ni ^y^a))- 



Notice that components of its position vector are the conserved densities appearing in Eqs. (1361) . 
(|37p . By construction, hypersurfaces Mp , Mp and M| have parallel tangent hyperplanes at 
the points corresponding to the same values of the parameters v l . Let us now apply an arbitrary 
affine transformation in A n+1 , x = Ax, where A is a constant (n + 1) x (n + 1) matrix. This will 
transform Eq. (13811 to 

x n+1 - x k v k + F Q (v) = (39) 

where the transformation v — > v will automatically be projective: 

_, iHv) - F n 



here l l , I are linear (inhomogeneous) forms in v. Since affine transformations preserve the 
properties of being parallel/tangential, they naturally act on the class of Godunov's systems. 

5.2 Godunov form of integrable quasilinear systems and generalized hyper- 
geometric functions 

Let H Slt ___ tSn+2 be the space of solutions of the system 

d 2 h Sj dh Si dh , ,„„, 

+ — •«-, hJ = l,...,n, i + j, (40) 



and 



duiduj Ui — Uj duj Uj — Ui du 



n+2 n 



d 2 h / y^ \ s i r s i v^ u j( u j ~ 1) dh 



^ Sj ) uAui - 1) uAui - 1) ^ 



duidui \ ^- ' J J Ui(ui-\) Ui{ui - 1) rr? u j ~ u % 9uj 

3 *■ .77^ 

v^ Sj Sj + s n+ i Sj + s n+2 \ dh 

^—f Ui — Uj Ui Ui — 1 / dui ' 



(41) 



for one unknown function h(u\, . . . , u n ). Here s\, ..., s n +2 are arbitrary constants. Elements of 
H S i,...,s n+ 2 are examples of the so-called generalised hypergeometric functions, see |18| . |29j . It 

is known that dim.ff aij .„ jgn+a = n + 1. Let go,gi, ■■■,gn be a basis of H Sl Sn+2 - Choose a time 

t q for each element g q of this basis. Here q runs from to n. It is known |27] . [29j that for each 
pairwise distinct q, r, s running from ton the system 

V^ / Uj(uj - l)u ijts - Ui(ui - l)u jjts , , , 

2^ {9q,Uj9r,Ui - fl^Sg.uJ ~^~ h O ■ [_g q g r ,uj ~ grg q ,u J )Uj,t a + 
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E .Uj(uj - l)u it t -Ui(ui - l)u jttq , , . ,, . 

{9r,Uj 9s,Ui - 9s, Uj 9r,u, ) h CT ■ \.9r9s,Uj ~ 9s9r, Ul )Uj )tq + (42) 
Uj — Ui 
l<j<n,j^j J 

E ,Uj(Uj-l)u it t r -Ui(Ui- l)u j)tr , , _ 

{9s,Uj9q, Ui ~ 5g,u,-5s,wJ h <7 ■ (ffa^g.Uj ~ 9q9s,uMj,t r = 0, 
Uj — Ui 

where j = 1, ..., n, <x = 1 + s% + ... + s n +2, possesses a disperssionless Lax representation and an 
infinity of hydrodynamic reductions. Moreover, any generic integrable 3-dimensional hydrody- 
namic type system with two unknowns is isomorphic to a system of the form (|42p . see |27| . It 
is also known that the system (l4"2|) possesses n + 1 conservation laws of hydrodynamic type, see 
|14j for n = 2 and |29| for general n. Therefore, it possesses a Godunov representation which 
can be constructed explicitly in the following way. 

Theorem 5 Let ho,hi, ...,h n be a basis o/i?2si,...,2s„+2- F° r eac h a ^ (3 = 0, l,...,n let f a p be 
a solution of the inhomogeneous linear system 

d 2 h 2si dh 2sn dh d^dvZ ~ &u~~5u~ . . -, . , . ,,„ N 

+ — --5— + , i,j = l,...,n, i^j, (43) 



duiduj U{ — Uj duj Uj — Ui dui Ui — u 

and 



d_h _( 1 + 2 Y^ \ 2si k\ 2si sr^ u j( u j- 1 ) dh 

duidui V ^-j" ° J Ui(ui — 1) Ui(ui — 1) *rf. Uj — Ui duj 



(44) 



y^ 2s j 2sj + 2s n+ i 2si + 2s n+2 \ dh 

^-^ Ui — Uj Ui Ui — 1 / dui 

dgadg^ _ dga:dgp_ dga dgp 

Edu, dUj d Uj dui Uj[Uj-l) s~d^9l3 ~ g^-9a 
1 TV - (- 1 + s l + ••• + s n+2) j 77 • 
. Ui - Uj Ui{Ui - 1) Ui(Ui - 1) 

Define new coordinates Vi, ..., v n and the functions F a p(v%, ■■■,v n ) by 

hi(ui,...,U n ) f a ,p (A r, 

ho(Ui,...,Un) n [ui,...,u n ) 

Then, in coordinates vi,...,v n , the system f^ffi) takes the Godunov form 



9F r A + (9F M \ | (MgA =Q i = 1 n (46) 

dv i )t q \ dv i)t r \ dv i)t s 

Proof: 

Substituting (|45p into (|46p and calculating derivatives of f a> p by virtue of (|43p . (j4"4"j) we obtain 
conservation laws of the system (|42p as found in |29| . More precisely, the left hand side of 
is equal to 



^Uj( U j-l)detW 3 

7 = 1 J J 
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where Rj is the left hand side of (|42p . the matrix W is defined as W = (w a ^) where wq 
hp(u\, ...,u n ), w a fi = Q-@- if a 7^ and W%j is the n x n-minor of W obtained by eliminatir 
column with hi and a row with -^—. In other words, the left hand side of (|46p is equal to 



-'.; 



so that (|4"6|) holds identically modulo (l42j) . 

Remark 1. The system (j46|) possesses a natural action of the group GL n+ \ x GL n+ \. Namely, 
the first copy of GL n+ \ acts on the times to, ...,t n , and in the same way on the basis go, ...,g n . 
This action corresponds to a change of basis go, ...,g n in H 81 Sn+a , The second copy of GL n+ i 
acts according to the Proposition 1. This action corresponds to a change of basis ho,...,h n in 

#2si,...,2s„ +2 - 

Remark 2. Note that F a g is defined up to an arbitrary linear combination of l,Vi, ...,v n and, 
therefore, f a R is defined up to an arbitrary linear combination oiho,hi,...,h n . This explains why 
f aj p satisfies a linear non-homogeneous system with the same homogeneous part as for elements 

from H2 Sl 2s n+2 - Moreover, the structure of the non-homogeneous part of the system (|43p . (j44j) 

containing linear combinations of -jSfr-fSf: — af^T^ anc ^ 'W^9p ~ ~5u~9 a ls dictated by the action 
of GL n+ i on the times to, ...,t n , and in the same way on go,...,g n , compare with (l4~2"jh 
Remark 3. It was proven in |27j that any 2-component integrable system is isomorphic to a 
member of the family (|42p with n = 2, or its appropriate limit. Therefore, Theorem 5 gives, in 
particular, a description of the Godunov form for any generic 2-components integrable system. 
If n > 2 there exist n-component integrable systems which do not belong to the family (|42p or 
its degenerations (see, for example, [30]). However, n-component systems constructed in [30J) 
have n conservation laws only and, therefore, do not possess a Godunov form. Probably, the only 
integrable systems possessing Godunov's form should belong to the family (|42p or its appropriate 
degenerations. 

5.3 Application to integrable potentials H(V, W) 

Let us now apply these results to the system (|lip . Set n = 2, choose a triple of indices s = 
0, q = 1, r = 2 and set to = — t, t\ = x, ti = y where t, x, y are the independent variables in 
pip . Note that the bases in iJ Slj ... jSn+2 and i?2si,...,2s„+2 are no longer independent since we have 
constraints on the functions Fi^, i*2,o> ^b,l; namely, F20 = V\F\2j ^0,1 = V2F1 2- To obtain (|lip 
it remains to set F12 = H, v\ = V, V2 = W. Therefore, /1 2 = Hho, /2,o = Hh\, /o,i = Hh>2- 
Moreover, we must have 

hi = a(u 1 ,u 2 )(g j gk,u 1 - 9k9j,ui) + b(ui,U2)(gjgk,u 2 ~ 9k9j,u 2 ) + c(u 1 ,u 2 )(gj, U2 9k,u 1 - gj, Ul gk,u 2 ) 

where i,j, A; is a cyclic permutation of 0, 1, 2. Indeed, hi must have the same structure in go> 9ii 92 
as coefficients at Ujti in (j42p . Substituting the expressions for ho,h\,h2 and /1 2, /2,o> /o,i hrto 
the equations for i?2si,...,2s„+2 an d (|43|l . (|44P . respectively, and using the equations (140 p . (|4ip for 
50)5'i)5 , 2) we obtain that si = S2 = —S3 = S4 = — i, along with the following expressions for the 
functions /io, /ii, ^2' 

ho{ui,u 2 ) = C[ — 2 —{gig2, Ul - 929i, Ul )+ 
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j(gi92,u2 - 929i,u 2 ) + 3(«i - u 2 )(gi,u 2 92,u 1 - 9i, ui 92,u2. 



hi 

hi(ui,u 2 ) = C( 

^u 2 



hi(ui,u 2 ) = C[ — ^—(g 2 go, Ul - go92, ul )+ 



(g2go, U2 - gog2, U2 ) + 3(«i - u 2 )(g 2 ,u 2 go, ul - g2, ul go, U2 ) 



U\ — 1 

h 2 (u 1 ,u 2 ) = C[ — ^-r(gogi, ul - 9i9o, Ul )+ 

\U 2 — L 

t(j9o9i,u2 - 9igo,u 2 ) + 3(ui - u 2 )(go,u 2 gi,ui - go, Ul gi,u 2 

U\ — 1 



where 



C=( Ul -lf/Z( U2 -lf/Z(u 1 -U 2 )- 1 /3. 

We also obtain the following system for H, 

U 2 -l "I" 1 

Ul c( Ul - U2 y U2 c{ U2 - Ul y 

The solution of this system reads 

'Ul — 1 



H(m,u 2 ) = G 

V u 2 - 1 - 

where G'(£) = 2 y 3 , , 2/3 , Summarizing, we obtain the proof of the following 

Theorem 4 Let the potential H(V, W) be defined parametrically in the form 

V u 2 - 1 J h ho 

where G'{t) = , 2 / 3( , ^2/3 and ho,hi,h 2 are three linearly independent solution of the hypergeo- 
metric system 

ui(l - ui)h UUUl - -uih Ul - -h = - h Ul + h U2 , 

6 9 6 u\ — u 2 6 u 2 — u\ 

-u 2 (l - u 2 )h U2:U2 - -u 2 h U2 - -h = h Ul + h U2 , 

S 9 d 111-M2 3 u 2 — ui 

7 ^ "-U2 "«1 

""111,112 TT 1 

6 u 2 — ui 

(this system coincides with ( (^0| j, fl^l| j /or £/je values of constants s\ = s 2 = —S3 = 54 = — | ). 
Then H(V, W) provides the generic solution to the system 1127]) . In particular, it does not possess 
any continuous symmetry from the equivalence group. 
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6 Appendix. The method of hydrodynamic reductions. Deriva- 
tion of the integrability conditions 

Applied to a (2 + l)-dimensional system ([1]), the method of hydrodynamic reductions consists of 
seeking multi-phase solutions in the form 

u(x,y,t)=u(R\...,R n ) 

where the phases R l (x,y,t) are required to satisfy a pair of (1 + l)-dimensional systems of 
hydrodynamic type, 

R\ = X l (R)R x , R l y = n l {R)R % x . 

Solutions of this form, known as 'non-linear interactions of n planar simple waves' [6b\ 135] , 
have been extensively discussed in gas dynamics; later, they reappeared in the context of the 
dispersionless KP hierarchy, see |19| |2"U] and references therein. Effectively, one decouples the 
(2+l)-dimensional system ([TJ into a pair of commuting n-component (l + l)-dimensional systems, 
which are called its hydrodynamic reductions. Substituting the ansatz u(i? , ...,R n ) into ([T]) one 
obtains 

(A*J„ - A - fi*B) diU = 0, i = l,...,n, (47) 

di = d/dR 1 , implying that both characteristic speeds X 1 and fi % satisfy the dispersion relation 

det(AJ„ -A-fiB) = 0, (48) 

which defines an algebraic curve of degree n on the (A, u)-plane. Moreover, A* and \x l are required 
to satisfy the commutativity conditions 



dj\ l djjji 



A^ - X 1 jm> - fi l 



(49) 



i ^ j, see |38] , In was observed in [13J that the requirement of the existence of 'sufficiently 
many' hydrodynamic reductions imposes strong restrictions on the system ([I]), and provides an 
efficient classification criterion. To be precise, we will call a system ([1]) integrable if, for any n, 
it possesses infinitely many n-component hydrodynamic reductions parametrised by n arbitrary 
functions of a single variable. Thus, integrable systems are required to possess an infinity of 
n-phase solutions which can be viewed as natural dispersionless analogues of algebro-geometric 
solutions of soliton equations. 

In this Appendix we illustrate the above procedure by characterising all Hamiltonians h(v,w) 
for which the system (J7|) is integrable by the method of hydrodynamic reductions. Rewriting 
Eqs. © as 

we seek n-phase solutions in the form v = v(R 1 ,R 2 , . . . , R n ), w = w(R 1 ,R 2 , . . . , R n ) where the 
phases (Riemann invariants) R l satisfy the equations 



R\ — X l (R)R l x , R y — tf{R)R x . 
The substitution of this ansatz into (|50p implies the relations 
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(51) 



here Vi = d{V, Wi = dtw, di = d/dR l . The condition of their non-trivial solvability implies the 
dispersion relation for A* and //', 

_\* 2 _i_/, \ i ,, i ^- r >h V- 4- h h — h 2 — n 

In what follows we assume that the dispersion relation defines an irreducible conic, which is 
equivalent to the requirement h ww ^ 0, h vv h ww — h% w ^ 0. Setting Vi = <ft l Wi we can rewrite 
(|5ip in the from 

6 % X = h vw 6 l + h. 



vwY ~ ""WW: 

i 



(52) 
A* = h vv 4> % + h vw n % 4> % + h ww [i l + h v 



b VW . 



We also require the compatibility of the relations x>i = (j^Wi, which gives 



b3 - 4> 
Expressing A* , fx 1 in terms of <j) % from ([52]) , 



Bdf dip 

didjW = — -A — —diW + — -A — —rdjW. (53) 



■2 

» i f^vw^P ~r "ww i "'WW "'vvY 

A — Ti ) M 



4> l (h vw (j) 1 + h ww y 

and substituting these expressions into the commutativity conditions (|49|) we obtain relations of 
the form dj(f) 1 = (• • • )djW, didjW = (. . . )diwdjW where dots denote certain rational expressions 
in 0*,0 J whose coefficients depend on the Hamiltonian density h(v,w) and its derivatives up to 
the order three. The compatibility conditions djdk<p l = dkdj<j) 1 take the form PdjwdkW = 
where P is a polynomial in fa, (j)j,(j)k- Setting all coefficients of this polynomial equal to zero one 
obtains the integrability conditions (Tl6|) . 
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